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Bouncing models have been proposed by many authors as a completion, or even as an alternative 
to inflation for the description of the very early and dense Universe. However, most bouncing models 
contain a contracting phase from a very large and rarefied state, where dark energy might have had 
an important role as it has today in accelerating our large Universe. In that case, its presence 
can modify the initial conditions and evolution of cosmological perturbations, changing the known 
results already obtained in the literature concerning their amplitude and spectrum. In this paper, 
we assume the simplest and most appealing candidate for dark energy, the cosmological constant, 
and evaluate its influence on the evolution of cosmological perturbations during the contracting 
phase of a bouncing model, which also contains a scalar field with a potential allowing background 
solutions with pressure and energy density satisfying p — we, w being a constant. An initial 
adiabatic vacuum state can be set at the end of domination by the cosmological constant, and an 
almost scale invariant spectrum of perturbations is obtained for ui ~ 0, which is the usual result 
for bouncing models. However, the presence of the cosmological constant induces oscillations and a 
running towards a tiny red-tilted spectrum for long wavelength perturbations. 

PACS numbers: 98.80.Cq, 04.60. Ds 



I. INTRODUCTION 

Bouncing models [ll, ^ have been widely investigated 
as a solution of the singularity problem, and possibly as 
an alternative to inflation as long as it can also solve, by 
its own way, the horizon and flatness problems, and jus- 
tify the power spectrum of primordial cosmological per- 
turbations inferred by observations. 

In the case where the contracting phase of a regular 
bouncing model is dominated by some matter content 
with a constant ratio between pressure and energy den- 
sity, p/e — w ~ const., it was shown by many authors, in 
different frameworks [3-8] , that this matter content must 
be dust-like, perhaps connected to cold dark matter, in 
order to obtain a scale invariant spectrum of scalar and 
tensor cosmological perturbations. 

On the other hand, since 1998 |9j], cosmologists were 
confronted with a highly unexpected observation: the 
Universe is presently in a state of accelerated expansion. 
This may be caused by the existence of some field violat- 
ing the strong energy condition, called dark energy, by a 
modification of general relativity at large scales, by the 
influence of some large scale inhomogeneities, or simply 
by a well suited cosmological constant. This last option 
is by far the simplest explanation to the present accel- 
eration of the Universe, although it poses a problem to 
quantum field theory on how to accommodate its ob- 
served value with vacuum energy calculations. Anyway, 
the so called ACDM standard model assumes that there 
exists a cosmological constant term in Einstein's equa- 
tions, which becomes dynamically important when the 
typical scale of the Universe has the size of the present 
Hubble radius. 

In bouncing models without a cosmological constant, 
vacuum initial conditions for quantum cosmological per- 
turbations are set in the far past of the contracting phase. 



when the Universe was very big and almost fiat, justify- 
ing the choice of an adiabatic Minkowski vacuum in that 
phase. However, if a cosmological constant is present, 
the asymptotic past of bouncing models will approach 
de Sitter rather than Minkowski spacetime. Further- 
more, the large wavelengths today become comparable 
with the Hubble radius in the contracting phase when 
the Universe was still slightly infiuenced by the cosmo- 
logical constant. Hence, the existence of a cosmological 
constant can modify the spectrum and amplitude of cos- 
mological perturbations. Note that this is not a question 
for inflation because initial conditions for quantum per- 
turbations and the moment of Hubble radius crossing in 
such models take place when the cosmological constant 
is completely irrelevant: the Universe is fully dominated 
by the inflaton field. 

The aim of this paper is to investigate this issue in de- 
tail in the context of a Friedmann-Lemaitre-Robertson- 
Walker geometry with a cosmological constant, and a 
scalar field with potential allowing a constant equation of 
state p = we for the background field, like the exponential 
potential in the scenario without cosmological constant. 
Hence, this paper can be considered as an extension of 
Ref. [1] through the introduction of a cosmological con- 
stant in the model. Here, as in Ref. our background 
scenario is not intended to be a fully realistic description 
of the contracting phase of a bouncing model, but to yield 
a suitable framework to calculate the spectrum of linear 
cosmological perturbations in bouncing models, and to 
study how it depends on the presence of a cosmological 
constant and on the equation of state of the matter con- 
tent. In our model, the bounce itself takes place at very 
short length scales, where the cosmological constant has 
no role. Hence, its presence does not modify the evo- 
lution of the background and perturbations in that pe- 
riod, and the descriptions provided in Refs. 0-0, liol - fl^ 
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can still be considered to be valid at the bounce. The 
main difference is originated from processes much before 
the bounce, when the initial conditions are set and the 
cosmological constant is not irrelevant. In that case, a 
Minkowski adiabatic vacuum can only be defined in a 
precise time domain, i.e., at the end of cosmological con- 
stant domination, but when the Universe was still very 
big and rarefied. However, even in this time domain, 
as the length scale associated with the cosmological con- 
stant, given by the present acceleration of the Universe, 
is not much bigger than the long wavelengths of physical 
interest today, the spectrum of these scales can still be 
slightly affected by the cosmological constant. And in- 
deed we will show, analytically and numerically, that the 
usual result for bouncing models, namely, that the fiuid 
should satisfy w ^ m order to have an almost scale in- 
variant spectrum of long wavelength perturbations, still 
holds, but the presence of the cosmological constant in- 
duces small oscillations and a small running towards a 
red-tilted spectrum for these scales. 

In the next section, we will present the background 
model and obtain the equations for the evolution of cos- 
mological perturbations on this background. In section 
III, we will discuss the choice of the initial state of the 
cosmological perturbations on this background. In sec- 
tion IV, we will obtain analytically and numerically the 
power spectrum of perturbations for the model presented 
in section III, and discuss its physical consequences. We 
end up with the conclusions. 



II. THE BACKGROUND MODEL AND THE 
EQUATIONS FOR SCALAR PERTURBATIONS 

The gravitational action we shall begin with is that of 
General Relativity with a cosmological constant, i.e. 




2A)d''a;, 



(1) 



where l^^ = (87rG„/3)^/^ is the Planck length in natural 
units {h — c — 1), and A is the cosmological constant. 

The geometry of the background is given by the spa- 
tially fiat homogeneous and isotropic line element in con- 
formal time: 



(3) 



where the potential energy density of the scalar field is 
given by 



U{ip) = Uosmh' (I), 



^ 3(1 - w)H^nA ^ I 1 
° 16ttG ' y6TrG{l + w)' 



(4) 



(5) 



and w is a constant. This potential was already studied 
in Refs. [Ulil for the case w = -1/3. 

In the case of a homogeneous and isotropic back- 
ground, one can find the scalar field solution 



6ttG{1 + w' 



■In 



tanh 



3(1 



(6) 

and under these conditions the energy density and pres- 
sure of the scalar field 



satisfy p = we (a prime denotes derivative with respect 
to conformal time). 

The Friedmann equations in conformal time read 



■H'-H^ = -iTrGa^{e+p), 



(8) 



(9) 



and e satisfies the conservation equation 

e' = -3n{e + p), (10) 



ds^ = a'^{7]){dr]'^ - S.jdx'dx^). 



(2) 



The matter content of the model is described by a 
canonical miminimally coupled scalar field ip with La- 
grangian 



where T-L = a' /a. 

In the present situation, where the pressure and energy 
density of the matter content satisfy p = we, with w 
constant, the solution for the scale factor in terms of 
cosmic time dt = a drj reads, 



a{t)^ao( — 



l/3(l+w) 



sinh — 



3^01(1 +a;)JTo^ 



2/3(l+w) 



(11) 
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where Hq = 72 km Mpc~^ is the present Hubble 
parameter, fJoi^ = eo/fcrit with ecrit = 3iJg/(87rG'), and 
Ha = a/ Hq . The subscript indicates the present vahies 
of the respective quantities. 

Note from Eqs. (|4l7p and p — we that ip — 
y/2{l + w)Uq/(1 — w) smhif/F, and the kinetic energy 
of the scalar field grows exponentially with Lp, as usu- 
ally expected for a scalar field in a contracting phase of 
a Friedmann model. However, the potential increases in 
the same way, and that is why p — we \s maintained. 
Nevertheless, there is the question about the instability 
against initial conditions of this tracking between poten- 
tial and kinetic energies of the scalar field in order to keep 
p = we m. the contracting phase, which is also an issue 
for Ref. Q, but note that it is not necessary that this 
tracking must be valid all the way to the bounce: here, 
as in Ref. Q , we are interested in the spectrum of long 
wavelength perturbations, which become comparable to 
the Hubble radius in the contracting phase at scales not 
much smaller than the Hubble radius today. Hence, our 
calculations for the spectrum of long wavelength pertur- 
bations should hold even if the tracking p = we ceases 
to be valid at smaller scales. This may affect the spec- 
trum of small wavelength perturbations, but it will not 
affect our forthcoming calculations for long wavelength 
perturbations if the tracking p = we holds at scales of 
the order of the Hubble radius today, when the contract- 
ing Universe is still very big. Hence, in the same way 
that the conclusions of Ref. [3| proved to be valid for 
more elaborate bouncing models, we expect that the re- 
sults we will present later on will also be valid in more 
elaborate models containing a cosmological constant. 

Note also that as l/-\/A is of the order of the Hub- 
ble radius today, then all scales of physical interest had 
causal contact in the contracting phase of the model be- 
cause the particle horizon in that phase is of the order of 
1/VvV. 

The evolution of linear scalar perturbations are de- 
scribed by the metric 



back I J 



(12) 



where represents the homogeneous and isotropic 

cosmological background given in Eq. ([2]), and the per- 
turbations hfj^y can be decomposed into 



hij = 2a'^{ip^^j ~ E^ij). 



(13) 



Due to the constraint equations present in the Ein- 
stein's equations, the evolution of quantum perturbations 
in a classical background is described by a single quan- 
tum field, the gauge invariant Mukhanov-Sasaki variable 
defined by (see Ref. ^15)] for details) 



a ( ^(^+ ^back^ 



(14) 



where Sep is the perturbed scalar field, and (fhack is its 
background solution. 

The Mukhanov-Sasaki variable satisfies the equation. 



V — V f V = U, 



(15) 



where 



o^y^inGNie+p) 
H ' 



(16) 



The equations above are not altered by the presence of 
a cosmological constant, essentially because Eqs. (I9ll0p 
are not modified by its presence and because, of course, 
6A = 0. 

In our choice of units a is dimensionless, hence we 
will define the dimensionless conformal time r) = ij/Rh, 
where Rh = l/(ao-ffo) is the co-moving Hubble radius. 
From now on we will omit the tilda over r/. We will 
also work with the dimensionless comoving wavenumber 
k = Rh/^, where A is the comoving wavelength of the 
perturbation modes. The region corresponding to long 
wavelengths today is the interval 1 < fc < 10'^. 

Taking the model with scalar field and scale factor 
given by Eqs. ^ and pT|) . respectively, one obtains that 



z(<) 



3(1 + w) a{t) 
2 cosh(7t) 



(17) 



and 



z a-Q 



(1 - 3w) 



sinh (7i) 



(1 -(- 3w) 



9(1 + ^)2 
2cosh^(7i) 



(18) 
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where 
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3^/^II(l + a;)iJo 



(19) 



For t —i' 0~ , 01 rj 0~ , one obtains from Eq. PT|) that 
a{t) (X t^/lscn-u))] Qj.^ conformal time, a{r]) cx 77^/(^+'^"''. 
This is the usual Friedmann evolution for p = we without 
a cosmological constant. In this regime, we obtain 



and a is given by Eq. ([TTjl. 

Solution pT|) is defined in two domains: — oo < t < 
and < t < oo. The first one describes a universe 
contracting from an asymptotic de Sitter spacetime in 
the far past to a singularity at t = 0. The second one 
describes a universe expanding from a singularity at i = 
to an asymptotic de Sitter expansion in the far future. 
Around t — 0, the field dominates the dynamics, and the 
cosmological constant is unimportant. These behaviours 
can be viewed by taking the limits e.g. , in the contracting 
solution, t —oo and i — > 0~ in Eq. (|lip . 

For t -oo, Eq. ([TTl) yields a{t) « exp(-\/Ai). In 
conformal time 



V + Voo 



_exp(^^^:^) 



(l+3w)/3(l+w) ' 



(20) 



where — t^oo < ?7 ^ 0, and ijoc is a positive constant, the 
scale factor behaves as 



0(77) 



a-o 



'^iv + Voo) 



(21) 



This is the usual de Sitter behaviour. In this case, the 
potential reads 



V{r,) = - « 

■z 4(?7 + r7oo)^' 



yielding the equation 



2 (9^' - 1) 



4(7? + ?7oo)^ 



Vk = 0. 



(22) 



(23) 



This equation is completely equivalent to an equation 
for a massive scalar field in a de Sitter spacetime, with 
mass given by 



3\/A 



(24) 



Its general solution reads 



(25) 

where the Hi, ' are Hankel functions of first and second 
kind, v — 3w/2. As fc(?7 + ?7oo) ^ 1, we can write this 
solution as 



Vk - ci(fc)(r7+r;„,)(i+3»)/2+C2(fc)(r7+,?oo)(i-3»)/2. (26) 



z 



2(1 - 3u;) 
(TT3z«)V' 



(27) 



In this situation, z oc a (see Eq. (jlZp ) and z" j z = a" /a. 

Note that the potential V{t) diverges to ±00 at the 
infinity past for w > 1/3 and w < 1/3, respectively, and 
diverges to ±00 near the singularity at i cx 77 ~ for 
w < 1/3 and w > 1/3, respectively. Hence, it must cross 
zero in the middle of the line — 00 < i < 0. In Fig. [T] we 
present the behaviours of the potential V{t) for t < in 
the cases w < 1/3, w = 1/3 and w > 1/3. 

Our idea is that the singularity at t = separating 
the contracting and expanding solutions can be elimi- 
nated through some new physics which produces a regu- 
lar bounce connecting these two phases. As in the region 
around t — the cosmological constant is unimportant, 
one can evoke the bounce descriptions provided e.g. , in 
Refs. [sI-ItI. [l0l - [l^ . For instance, the quantum cosmologi- 
cal bounces with a perfect fluid studied in Refs. [tI. [l0l - [l3 | 
present a regular scale factor given by 



a(T) = flb 



(28) 



where dT 



,l-3uj 



dr], and Ub and Tf, are positive con- 



stants. Note that for |r| Tb, this solution approaches 
the classical Friedmann solution for a perfect fluid given 
by the limit t -J> 0" of Eq. ITT]): a{ri) cx . Hence, 

the scale factor (fTTj) can be smoothly connected to the 
scale factor (p8|) . 

It was shown in Ref. [l^ that the potential present in 
the equations for the perturbations around these quan- 
tum bounces reads 



ViT) ^ - = 



2 a 



3(l-«i) 



3(i-z.)r,^ 



-,l-|-3'!I> 



2 3(1-1 



1 



(29) 

where Eq. (|28l) has been used. Their shapes are presented 
in Fig.[5]for the cases ?i'<l/3, ?i' = l/3 and w > 1/3, and 
they tend to the potential V{t) of Eq. ^ for 7t << 1 
in the limit |r| ^ Tb, but they do not diverge in ry = 0. 
Hence, when these quantum effects become important 
near 77 = inducing the bounce, the two disjoint parts of 
the classical potentials presented in Fig. [T] corresponding 
to the contracting and expanding classical universes sep- 
arated by a singularity can now be softly connected with 
the potentials presented in Fig. [5] Then one can evolve 
smoothly the perturbations from the contracting phase 
to the expanding phase, and calculate their properties in 
the present era. For other regular bouncing models, the 
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FIG. 1: Behaviour of the potential V{t) given by Eq. (fTS)) for FIG. 2: Behaviour of the quantum bouncing potential given 
three different values of w. by Eq. ([28|l for three different values of w. 



situation must be similar, although somewhat more in- 
trincate in the case there is an extra field which induces 
the bounce. 

In order to accomplish this program, one must set the 
initial conditions for the perturbations in the past con- 
tracting phase. Without the cosmological constant, the 
Universe tends, in the far past, to Minkowski spacetime, 
where the potentials become null (see Fig. [2]). Hence, 
an adiabatic Minkowski vacuum can be prescribed there. 
In the presence of the cosmological constant, neither the 
Universe tends to Minkowski spacetime in the far past 
(in fact, it tends to de Sitter spacetime), nor the poten- 
tial becomes null there (see Fig. (TJ except for w = 1/3, 
which is not physically interesting because it does not 



yield an almost scale invariant spectrum of perturba- 
tions). However, as we have shown above, the poten- 
tial crosses zero somewhere in the middle of its evolution 
(which coincides with the moment when physically in- 
teresting long wavelengths perturbation modes become 
smaller than the Hubble radius), and perhaps one could 
define an adiabatic Minkowski vacuum there. We will 
show in the next section that this is indeed possible for 
the scalar field model presented above. 
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III. THE CHOICE OF INITIAL STATE FOR 
THE QUANTUM PERTURBATIONS 

In this section we will check whether an adiabatic 
Minkowski vacuum can indeed be prescribed in the time 
interval when the potential ((TH)) is negligible and the 
Universe is starting to be dominated by the field. In 
this regime the scale factor is approaching the form 
a{r]) cx 7^2/(1+3^)^ Taking flo^ w 0.3 and fl^ « 0.7, 
the zero of this potential occurs when y = yy = "fty 
(see Eq. (IT^ for the definition of 7). As yy depends on 
the w parameter, we will consider the following values for 
uj < 1/3: 



Vk(x) 



1 



■ exp 



where flk (x) must satisfy the equation 



2flk dx^ 



and fl = k'^+ (3x. 

Order by order, one has: 



3 / rfilfe 
dx 



(34) 



77^ ^ ' (35) 



yy(w « 0) 
yy(w« 1/8) 
yvii^ ~ 1/4) 



-0.34 , 
-0.23 , 
-0.13. 



Around these points, one can numerically approximate 
Eq. (|15p . now expressed in terms of the modes Vk{r/), 
through 



dx^ 



+ (fc2 + I3x)vk = 0, 



where 



yielding 



3(1 w)vTII / a dV\ 
2 \ao dy J 



(30) 



(31) 



Piiu « 0) 
l3{oj « 1/8) 
(3iu; « 1/4) 



-1.05 
-1.65 
-2.31. 



and X = — rjy, with rjv being the dimensionless confor- 
mal time corresponding to yv defined above. 

It is important to remark the dependence of /3 with 
17 A by looking at Eqs. PTI) and ([TH]). Noting that yv is 
independent of J^a, and as flaw + = 1, one gets 



(32) 



where /3o.7 are the values of /3 for n\ — 0.7. One can see 
that /3 ->■ as ilA ^ 0. 

Note also that, although yv is independent of f^A, Vv 
depends on fl^ as 



qj(1+3w)/[6{1+w)]q^^1/[3{1+w)] 

i/m+w)] W(o.7), (33) 



(l+3ii.)/[6(l+«))] 



where, again, r]v{o.7) are the values of rjv for 0,^ — 0.7. 
One can see now that ryy — >■ 00 as ^a 0, as expected. 
In this last calculation, we have assumed that the field 
dominates at rjv- 

The adiabatic vacuum is defined by the solution 



(0)^2 



16 A' 



16/6 



490 
256 



(36) 



where the upper indices {n) denote the order of the ap- 
proximation. Hence, an adiabatic Minkowski vacuum 
can be obtained if the parameter expansion /3^//^ sat- 
isfies /3^//^ << 1. In fact, as (i is of order unity, 
x << 1, and the long wavelengths of physical interest sat- 
isfy 1 < fc < 10^ the condition /3V/fc ~ P'^/k^ « 1 is 
satisfied. Note that for the largest scales [k approaching 
1), deviations from the Minkowski vacuum become more 
significant, and one should expect modifications against 
the standard results. 

This problem can be presented under another point 
of view. A Minkowski vacuum can be defined for quan- 
tum perturbations with wavelengths much smaller than 
the Hubble radius, defined by TZnit) = 1/H{t). From 
Eq. (fTTj). one obtains that 



7^i^(^) 



-tanh(-j/). 



(37) 



One has to compare this quantity with the physical 
wavelength Aphys = a A which, from Eq. pT|) . reads 



\ \ now 

^phys — Aphys 



l/3(H-w) 



sinh 



2/3(l+Lj) 



i-y)- (38) 



The maximum value of TZh, at i — >■ — 00, is W]^^^{t) — 
A"-'^/^, while Aphys diverges there. Eventually, they can 
be comparable at some time in the contracting phase. 
As in the case this is true one expects to obtain a similar 
spectrum as the one obtained in bouncing models with- 
out a cosmological constant, we will concentrate on the 
case < w << 1, which yields an almost scale invariant 
spectrum. The quantities defined in Eqs. (|37l) and 
are comparable when 



^1/3^1/6 V^yi 

Uoj A 7-7 now 



sinh"'^^'^(— y) 
cosh(-?/) 



(39) 
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where 



is the Hubble radius 



sinhi/3(-y)/cosh(-y) < 0.73, f^o™^^^' 



today. As 
'OA ~ 0.63, and 

10-3 < A"^^p;/7^^°* < 1, this equality can be achieved 
for some finite domain of y. 

Note also from Eq. p9l) that this domain interval of y 
could be extended to large values oi\y\ if f^A were much 
smaller than our prescribed values. This can also be seen 
from the analysis coming from the potential, where a 
smaller f^A would result in a smaller /3 in Eq. and 
an adiabatic vacuum could also be achieved for smaller 
values of k. 



Hence, substituting the zero order term given in 

Eq. (ISni) in the solution dM]) of Eq. ((201), we obtain 



Vk{x) 



1 



■ exp 



2(fc2 + /3a;)i/4 
and the initial conditions are given by 
1 



«fc(0) 



2^/fc 



exp 



(40) 



(41) 



IV. SPECTRUM OF QUANTUM 
COSMOLOGICAL PERTURBATIONS 

Let us now calculate the spectrum of quantum cos- 
mological perturbations for this scenario. In Section HI, 
we have shown that an adiabatic Minkowski vacuum, for 
the case of a canonical scalar field, could be prescribed in 
the time interval where the potential becomes negligible 
and the Universe is starting to be dominated by the field. 



J 



dx 



«-.(0)(^+.fc). 
We have, therefore, to solve 



(42) 



v'l+\k' Ufc = 



with initial conditions given by Eqs. and 



(43) 



We calculated numerically the solution of Eq. (|43p by changing the time variable from 77 to y, defining a new 
function Uk = a^^'^Vk, and setting the above initial conditions at yini = yv- Taking the potential (jlSp . the transformed 
equation reads 



dy^ 



4fc2 



w 



w 



9(l + u;)2[l](^i+3"')/2no^]2/[3(i+-)]sinh4/[3(i+'")](_y) {I + w)^ {I + w)^ smh\y) cosh2(y) 



The results are shown in Fig. [3l 

The solutions of equation can be expanded in powers of k^ according to the formal solution [If 



Ai(fc) 



1-e 



dri 
(77) 



(77) d77 ' 



-42 (fc) 



d77 
(77) 



drj 
(77) 



Uk = 0. 

(44) 



, (45) 



In Eq. (j45|) . the coefficients ^i(fc) and A2{k) are two 
constants depending only on the wavenumber k through 
the initial conditions. 

Since expansion (j45p is valid at all times during con- 
traction, the Ai (k) and A2 (fc) dependencies coming from 
the initial conditions hold when the Universe enters the 
field dominated phase just before performing the bounce. 
From then on, the evolution is guided by the particu- 
lar physics of the bounce. For instance, in the quantum 
bounce with potentials shown in Fig. [21 everything goes 
as described in Ref. 0, with power spectrum in the ex- 
panding phase given by 



V (X k^\A2{k)\^' 



(46) 



and spectral index 



dlnjV) 
dln(fc) ■ 



(47) 



This is a general feature of bouncing models: the spec- 
trum A2{k) of the growing mode of the contracting phase, 
which will be the decaying mode of the expanding phase 
after the bounce, is transferred to the growing mode of 
the expanding phase and dominates over Ai(k). Hence, 
to find the spectrum after the bounce one has to obtain 
A2{k), which, in the case without cosmological constant, 

3(^-1) f-i 

reads A2{k) cx fcso^+D (see e.g., Ref. (7|), yielding the 
spectral index 



12a; 
l + 3a;' 



(48) 
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FIG. 3: Numerical results of Eq. (gU for w = 10"^, 1/8 and 
1/4. Each curve shows Uk{y) as a function of y for different 
values of k, as indicated in the graphs. 



— l\ h KXf 

V 3/3 



(49) 



while its first derivative reads 



dvk 

dx 



X exp 



3/3 



(50) 



The second approximate solution comes from the re- 
mark that at rj^ and afterwards, up to the bounce phase, 
the evolution of the background is dominated by the 
scalar field, where the potential approaches the form 
given in Eq. (I27[) . yielding the solution 



Vkiv) = ^ [Ci{k)Hlikr^) + C2{k)Hl{k7^)\ , (51) 

where v = 3(1 — w)/2(l + 3u'). In the domain |fc?y,| >> 1 
this solution reads 



and 

dvk 

dx »)=77, 

where 



Biik) 



2Vk 



2Vk 



2ikrj^ 



B2{k) 



2Vk 



2ikr]t 



(52) 



{Bi{k) 



ik - 



2??* 



B2{k) 



ik - 



277* 



Bi{k) = 2Ci(fc)exp i^A;?7» 
B2{k) = 2C2(fc)exp[-z(fc?7*-^- J 



TTl' TT 
T ^ 4 



ai 



T{v + 3/2) _ 2{l-3w) 
r{iy - 1/2) ^ (1 + 3w)2' 



(53) 

(54) 
(55) 

(56) 



Performing the matching between Eqs. (I49l50p and 
Eqs. (I52I53P at 77* , one gets for Ci(fc) and C2(fc), 



Ci{k) 



exp 



/2k^ 
-ii^—+kx^+a2 



(57) 



with w ~ for an almost scale invariant spectrum. 

In order to analytically predict the behaviour of the co- 
efficient A2{k) for the case with a cosmological constant, 
we have taken two approximate solutions of Eq. (l43t close 
to the point where the potential vanishes, where it has 
the form given in Eq. pop . and matched them at a point 
X* = (?7, — r]v), where < << 1. 

The first approximate solution comes from Eq. (1401) . 
and is given by 



and 



C2(fc) 



1 



1 



1L\ 



2 

X exp 



./2fc3 



where 



a2 = kf]^ 



(58) 



(59) 



In the limit krj , just before the new physics which 
generates the bounce, the solution Vkirj) approximately 
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reads 

' krj\'^ 1 




2 ) T{v + 1) 

+ i[Cx{k) - C2{k)]cot{vn) 

' kT^Y" i[C2{k) - C^{k)] 
r(l — v) sin(t^7r) 



Ci(k) + C2{k) + 



(60) 



and the coefficient A2 (fc) of the growing mode of the con- 
tracting phase is given by 



l/fc\-2>- f 

4 V2/ [r(f - v) sin(zy7r)]2 



l + _cos(2a2) + ^-^3 



Caculating the spectral index as defined in Eq. (|47l) . 
we find 



ns = 1 
2 



f2u; 



f + 3w 



k 



sin(2Q;2) 



^ + ^(f + 2a;,) cos(2a2) + ^ cos(2a2) 
4?7, 4 4 



sin(2a2) 



(62) 



Substituting the parameters it; 0, /3 w — f.05, |ry, | ~ 
|?7y| « 2.19, << 1, ai 2, a2 ~ fc?/* — tt, and 1 < 

< 10^ in Eq. we obtain an almost scale invariant 
spectrum. Besides the usual 12w/(l + 3w) result, there 
are additional terms in Eq. (j62l) inducing a running red- 
tilted spectrum and oscillations, both decreasing with k. 
Note that for a vanishing cosmological constant we have 
/3 « and [t;*! — >■ 00. In this case, the extra terms in 
Eq. (j62p disappear and — > 1 even for small values of 
k. 

In order to check numerically this analytic calculation, 
we took the following steps: from the numerical solutions 
Uk = a^^'^Vk presented in Fig. [31 we obtained Vk, evalu- 
ated it at very small y {y ~ — 10~^^), expressed the re- 
sult in conformal time (whose relation with y is trivial at 
field domination), multiplied the result by ri~'^~^/'^ (see 
Eq. (|60p ). and differentiated the final result with respect 
to 77 in order to isolate A2{k). The results are shown in 
Figs. 12 and El 

It can be seen that A2{k) and ns follow the predicted 
behaviour, A2{k) cx fc3(™-i)/2(i+3») and ns~l + 
for > 1. We have also verified that small oscillations 
and a red-tilted running with amplitudes decreasing with 
k are superimposed to the power law overall behaviour, 
as predicted in Eqs. (|6T1) and ([62l) . Note from Fig. 5 that 
the oscillations do indeed become smaller for smaller JIa , 
showing that they are a consequence of the presence of 
the cosmological constant. 



10 



0.1 



0.001 



10-- 



10- 



10 



50 100 



500 1000 



FIG. 4: Numerical results for the behaviour of |A2(fc)| for 
w = 10"^, 1/8 and 1/4 evaluated at j/ = -10"^^ The solid 
lines show the numerical results and the dotted ones show, 

3(-^-l) 



sin(2a2) -(61) for comparison, a curve proportional to fc^C 



3u+l) 



1.016 - 




FIG. 5: Numerical results for ns{k) evaluated sd, y — — 10~^^, 
obtained using w = 10"'^. The solid line indicates the result 
obtained using f^A = 0.7, the dashed line for f^A = 10""^ and 
the dotted line for f^A ~ 10~®. Note that the oscillations 
become smaller for smaller f^A, showing that they are due to 
the presence of the cosmological constant. This result is in 
agreement to Eq. H62[) . 



V. CONCLUSION 

In this paper we have investigated the effects of the 
presence of a cosmological constant in the contracting 
phase of a bouncing model. It turns out that the ini- 
tial vacuum state, usually determined in the contracting 
phase when the Universe was very large and rarefied, is 
affected by the presence of the cosmological constant. In 
order to get an almost scale invariant spectrum, one still 
must have some phase with dust-like contraction, but 
now the spectral index gets a red-tilted running and os- 
cillations directly caused by the cosmological constant. It 
is interesting to realize that bouncing models allow such 
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an important role to the cosmological constant in the 
physics of primordial cosmological perturbations, which 
is not at all the case for always expanding models. This 
opens a new area of research, which is to investigate the 
influence of other models of dark energy on the primor- 
dial spectrum of bouncing models.. In other words, if 
the Universe had really bounced in the past, investigat- 
ing its primordial spectrum can yield information about 
dark energy. 

There is also the question about the possibility of an 
enormous growth of perturbation amplitudes in the de- 
flationary contraction in the far past of the model, as 
discussed in other contexts [l3|. Note, however, that the 
cosmological constant in our model is small and this al- 
most de Sitter deflation will take place when the Universe 
was very large compared to the present Hubble radius, 
and for a fixed time interval. In conformal time, this time 
interval can be estimated using Eq. (I^H)) by saying that 
it should be smaller than 77, -I- 7700 given by this equation 
when t « — 1/a/A. Taking the usual values 51 a =0.7 and 
^Ow = 0.3, one can see that 77* -I- 7700 < 1- Note that 
for a cosmological constant dominated model {^ow ~ 0), 
one would obtain 77, -I- 7700 >> 1. Let us examine the 
behaviour of the Bardeen potential in this phase. From 
Eq. ([25)1 we obtain that 

for w = 0. Hence, as this almost de Sitter deflation- 
ary phase will not take long enough in conformal time, 
V*+Voc < 1, because of the smallness of the cosmological 
constant, perturbations will not grow alarmingly in this 
epoch. Now, once the Universe leaves this deflationary 
contraction to a non-deflationary contraction when it is 
still very large, then it can be subjected to dissipation 
effects, as the ones discussed in Ref. 0, which could dis- 
sipate the existing inhomogeneities. Another approach to 
this problem should be to think in terms of the Anthropic 
Principle and state that the Universe is composed funda- 
mentally by a small cosmological constant (dark energy) 
and some matter content as the one used in our model 
(dark matter?). In many regions it will expand to de 



Sitter and it will freeze, in some it will contract inhomo- 
geneously, and in a few homogeneous regions within one 
particle horizon size it may contract to make a bounce, 
where new particles (photons and baryons) will be cre- 
ated, and expand again to a Universe with some galaxies 
and stars where intelligent life can exist. The results of 
our paper can then be applied to this last possibility, the 
only one which can interest us. Of course these tenta- 
tive answers to this basic question must be worked out 
more precisely, but we think that a final and complete 
answer to the issue on why the primordial Universe was 
so homogeneous, in any cosmological scenario, demands 
a theory of initial conditions, perhaps quantum cosmol- 
ogy. One interesting investigation should be a quantum 
cosmological analysis of eternal asymptotically (in time) 
de Sitter models. 

Our approach here was concentrated on general fea- 
tures of the spectrum of cosmological perturbations in 
the presence of a cosmological constant in general bounc- 
ing models, and because of that we were not able to fix 
the amplitude of the perturbations. Our next step will 
be to take a specific model in which the physics at the 
bounce fixes the bounce scale, and hence the amplitude 
of the perturbations, in order to determine the influence 
of the new fetures of the spectrum of primordial pertur- 
bations we obtained in this paper in the anisotropics of 
the microwave background radiation, and to compare the 
results with observations. Another interesting problem 
should be to investigate the situation where the required 
dust-like contraction was not caused by a scalar field but 
by a hidrodynamical fluid with cj. = w. In this case, the 
prescription of an adiabatic vacuum can be much more 
involved because of the smallness of the sound horizon. 
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